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ABSTRACT 

We have analyzed the phase space distribution of a sample of about 900 non-kinematically selected low metallicity stars in the solar vicinity. 
The stars primarily represent the thick disk and halo populations of the Milky Way. We aim to identify overdensely populated regions in phase 
space, which we interpret as signatures of star streams passing close to the Sun. The search was conducted in a space constructed from the 
angular momenta and eccentricities of the stellar orbits. Besides recovering all well known star streams in the thick disk, we isolated four 
statistically significant phase space overdensities amongst halo stars. One of them is associated with a previously known halo star stream, but 
three of them are novel features, which we propose be also considered as genuine halo streams. 
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1. Introduction 

The recent large scale surveys of the stellar population of the 
Milky Way have revealed the presence of long coherent stel- 
lar streams. Most spectacular are the Sagittarius stream, which 
wraps around the entire Galaxy nearly perpendicular to the 
midplane, (Ibata et al. 2001, Martinez-Delgado et al. 2001, 
Majewski et al. 2003, Belokurov et al. 2006), and at low 
Galactic latitudes the Monoceros stream (Yanny et al. 2003, 
Ibata et al. 2003). These are interpreted as tidal debris associ- 
ated with disrupted satellite galaxies, which were accreted by 
the Milky Way. Numerical simulations of such merger events 
have shown that such debris streams survive as coherent struc- 
tures over gigayears so that they are witnesses of the formation 
history of the Galaxy (Helmi 2004, Law, Johnston & Majewski 
2005, Penarrubia, Marti'nez-Delgado & Rix 2005). 

Stellar streams have also been discovered as overdensities 
in the phase space distribution of stars in the solar neighbour- 
hood. Helmi et al. (1999) detected signs for a star stream in 
the halo population passing close to the Sun which was later 
confirmed by Chiba & Beers (2000). Prominent streams or 
more correctly 'moving groups' in the disk population of stars 
in the vicinity of the Sun are the Hercules, Pleiades-Hyades, 
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and Sirius streams (Eggen 1996, Dehnen 1998). These moving 
groups are very probably not related to tidal debris streams, but 
originate dynamically within the disk itself. Dehnen (2000) has 
demonstrated that stars in the Hercules group are presumably 
in outer Lindblad resonance with the inner bar of the Milky 
Way which leads to focussing of stellar orbits into certain parts 
of phase space. Quillen & Minchev (2005) have suggested that 
the other moving groups are due to resonances with spiral den- 
sity waves. All these features can be traced both in the thin and 
the thick disk populations. The Arcturus stream (Eggen 1996), 
on the other hand, is attributed to the thick disk only and is 
probably of extragalactic origin (Navarro, Helmi & Freeman 
2004). Arifyanto & Fuchs (2006) have analyzed the kinematics 
of nearby subdwarfs and identified another stream in the thick 
disk which resembles closely the Arcturus stream. The same 
feature has been detected independently by Helmi et al. (2006) 
in the phase space distribution of nearby F and G stars. The na- 
ture of this stream is at present unclear Helmi et al. (2006) ar- 
gue in favour of an extragalactic origin, because they find that 
its phase space structure resembles that of numerically simu- 
lated debris streams of disrupted satellite galaxies. On the other 
hand, it could as well owe its existence to similar dynamical 
effects which caused the Hercules group. If the members of 
the Hercules group are at outer Lindblad resonance with the 
Galactic bar (Dehnen 2000), the members of the new group 



2 



C. Dettbam et al.: Signatures of star streams in the phase space distribution of nearby halo stars 



Fig. 1. Slices through the phase space distribution of halo stars in the solar neighbourhood. The panels a, b, c, and d show the 
velocities of stars on orbital planes inclined relative to the direction towards the North Galactic Pole by angles of v = 0° - 45°, 
45° - 90°, 90° - 135°, and 135° - 180°, respectively (top left to bottom right). Vaz denotes the tangential velocity of the stars and 
Vae is a measure of the eccentricities of the stellar orbits (cf. section 2). 



would be at the outer 4: 1 ultraharmonic resonance with the bar 
(Fuchs, unpublished). These questions can only be settled, if 
the elemental abundance pattern of the members of the streams 
are known. They would indicate that either the stars belong to 
the field population of stars in the Milky Way or to the popula- 
tion of a disintegrated dwarf galaxy with different abundances 
(Williams & Freeman, in preparation). 

In this study we investigate the phase space distribution of 
halo stars in the solar neighbourhood, searching for signatures 
of streams in the Galactic halo. 

2. Search Strategy for Star Streams in the Halo 
Population 

Our analysis is based on the sample of non-kinematically se- 
lected low metalUcity stars of Beers et al. (2000). We have 
drawn from that catalogue the space velocity data of all stars 
with metallicties [Fe/H] < -1, which leads to a sample of 
897 stars. To the heliocentric space velocity vector {U, V, W) 
of each star, with U pointing towards the Galactic center, V 
into the direction of Galactic rotation, and W towards the North 
Galactic Pole, respectively, the velocity vector of the motion 
of the Sun relative to the local standard of rest {Uq, Vq, Wq) = 
(10.0, 5.3, 7.2)km/s was added (Dehnen & Binney 1998). Next, 
the velocities are transformed from the reference frame of the 
local standard of rest to an inertial frame by adding to the V ve- 
locity components the circular velocity of the local standard of 
rest around the Galactic center, for which we adopt the value 
of Vlsr = 220 km/s. In the following we assume a spherical 
Galactic gravitational potential and neglect the flattening of the 
disk potential and any asphericity of the dark halo potential. 
Chiba & Beers (2000) have modelled the Galactic gravitational 
potential with aspherical Stackel potentials. They showed that 
the distribution of halo stars in the phase space spanned by the 
isolating integrals of motion of the stellar orbits in such a poten- 
tial can be closely mapped into the corresponding phase space 
constructed assuming a spherical potential, which justifies our 
procedure. 

In a spherical potential a star moves in a constant orbital 
plane. Most of the stars which we have extracted from the 
Beers et al. (2000) sample have distances less than 1 kpc from 
the Sun, and for the purpose of calculating orbital parameters 
may be assumed to be at the Sun. The azimuthal velocity of 
a star is obviously given by Vaz = + ^lsr)2 + and 
the angular momentum is then L = Rq ■ Vaz, where de- 
notes the galactocentric distance of the Sun for which we as- 
sume /?0 = 8 kpc. The inclination angle of the orbital plane 
relative to the direction towards the North Galactic Pole is 
given by V = arctan((V -i- Vi^sr)/W). We restrict in our anal- 
ysis the inclination angles to the range 0° < v < 180°. Stars 
in an orbital plane with an inclination angle larger than 180° 
are treated as retrograde stars in the orbital plane with an in- 



clination angle v - 180°. The aim of our study is to identify 
star streams among the halo stars. Thus we are searching for 
groups of stars on essentially the same orbits which pass at 
present through the solar neigbourhood. For this purpose we 
have generalized the search technique of Arifyanto & Fuchs 
(2006) and have grouped the stars in our sample according to 
the inclination angles of their orbital planes. In each v-slice 
we have then searched for overdensities in a space spanned by 
angular momentum, which defines the mean guiding center ra- 
dius of a stellar orbit, and eccentricity, which describes the ra- 
dial excursion of the orbit from its mean guiding center radius. 
Since all stars are at present close to the Sun, they all have the 
same phase along their orbits. Using the theory of Galactic or- 
bits by Dekker (1976) Arifyanto & Fuchs (2006) have shown 
that the search can be conducted in practice in the space of 
Vaz and Vae = V^^^ + 2(Vaz - Vlsr)^. Vaz is directly related 
to the angular momentum and Vae, defined by the diff'erence 
between the actual orbital energy and the orbital energy of a 
star on the circular guiding center orbit, is a measure of the ec- 
centricity of the orbit, e - Vae/ V2yLSR- The latter assumes a 
flat Galactic rotation curve. Vae is also related to the radial ac- 
tion integral of the orbit and should be quite robust, even if the 
Galactic potential underwent changes in the past. In order to 
enhance the overdensities of stars in phase space we have per- 
formed a three-dimensional wavelet transform with a Mexican 
hat kernel •7C(x - x') = (3 - (x - x'f/a^) exp -((x - x'f/2 a^) 
(cf. Skuljan, Hearnshaw & Cottrell 1999) of the V^z - Vae - v 
data cube. After some experimentation we adopted a wavelet 
length scale of a = 20 km/s which leads to the clearest results. 
The wavelet transform conserves the total size of the occupied 
phase space, because the volume integral over TC is equal to 
zero. 

3. Results and Discussion 

3.1. Star streams in the halo population 

In Fig. 1 we show scatter plots of slices through the three- 
dimensional phase space distribution of the stars in our sample. 
Each sUce is 45° wide, and covers 0° to 45°, 45° to 90°, 90° 
to 135°, and 135° to 180°, respectively. Corresponding slices 
through the three-dimensional wavelet transforms of the data, 
which enhance and delineate the overdensities in a quantita- 
tive way, are shown in Fig. 2. A dififerent projection of the 
three-dimensional data cube onto inclination angle v versus 
azimuthal velocity Vaz is shown in Fig. 3. Since the V-shaped 
features in the Vaz - Vae planes are rather narrow, we have in- 
tegrated in Fig. 3 all phase space densities over Vae • A three- 
dimensional rendering of the data cube is shown in Fig. 4. As 
can be seen from Figs. 2, 3, and 4 there is a lot of fine structure 
in the phase space distribution of the stars. In the most fully 
populated part of phase space the inclination angles of the or- 
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Fig. 2. Slices through the wavelet transform of the phase space distribution of halo stars in the solar neighbourhood constructed 
analogous to the slices in Fig. 1. The phase space densities are colour coded at 10 % levels according to a Unear colour table 
ranging from red, over yellow, green, blue to dark violet. The labels indicate halo star streams. 



bital planes are about 90 degrees. These stars have disk-like 
kinematics and most of them are even at these comparatively 
low metallicities members of the thick disk of the Galaxy. The 
diagrams in the upper right and lower left panels of Fig. 2 show 
clearly the Hercules stream at Vaz = 170 km/s, the stream found 
by Arifyanto & Fuchs (2006) and Helmi et al. (2006) at V^^ = 
120 km/s, signs of the Arcturus stream at Vaz - 100 km/s , and 
possibly a further feature at Vaz = 75 km/s, which will be in- 
vestigated more closely in a separate study (Klement et al. , in 
preparation). 

However, our main interest here is star streams in the gen- 
uine halo star population. Tremaine (1993) has estimated a fill- 
ing factor by star streams in the halo around the Sun of about 
1, and we expect to find a few genuine halo star streams in the 
sample. In order to isolate these we have indicated in Fig. 3 the 
area populated predominantly by stars with disk-like kinemat- 
ics. The dashed Unes indicate the loci of stars with W = ±50 
km/s, respectively, and V^z ranging from 20 to 270 km/s, which 
should bracket the velocity range of thick disk stars. Outside 
this range we find at velocities Vaz > three distinct iso- 
lated features which we propose to interpret as signatures of 
star streams passing through the solar neighbourhood. Two of 
them (labelled HB and Si) are clearly visible on the right side 
of the bottom right panel of Fig. 2 and in the upper right part 
of Fig. 3. Their centroid positions in phase space are summa- 
rized in Table 1. Interestingly the feature labelled HB is the 
stream detected by Helmi et al.(1999). Chiba & Beers (2000) 
confirmed this discovery later using the same data of Beers et 
al. (2000) as analyzed in this study, but their 'clump stars' were 
isolated in a different way from our method. Actually we draw 
confidence in our search technique from the fact that it recov- 
ers a previously known star stream in the halo population. The 
stars in the third feature at Vaz = 140 km/s and v s; 30° have 
the same orbital parameters (cf. section 3.3) as the stars in the 
stream newly found by Arifyanto & Fuchs (2006) and Helmi 
et al. (2006). The members of this stream have disk-like kine- 
matics, but comparatively large vertical W velocities. Since the 
incUnation angles v depend on the vertical W velocity com- 
ponents, such stars will appear as linearly smeared features in 
Figs. 3 and 4. The other linear structures in Figs. 3 and 4 are, in 
our view, also due to this efl'ect. The features in the right lower 
half in Fig. 3 are therefore very unlikely to be associated with 
a genuine halo stream. 

Stars at negative azimuthal velocities are on retrograde or- 
bits and belong exclusively to the halo population. It is more 
difficult to discern overdensities in this part of phase space, 
because the V-shaped branches in the Vaz - Vae diagrams in 
Fig. 2 become very narrow and the inclinations of the orbital 
planes are evenly distributed over v in Fig. 3. However, we find 
two further isolated features in this part of phase space. One 
is located in the upper left frame of Fig. 2 at negative Vaz dis- 
tinctively above the narrow Vaz - Vae line of the bulk of the 
stars. This feature is also seen very clearly in the left side of 



Fig. 3. Projection of the wavelet transform of the phase space 
distribution of halo stars onto a plane spanned by the inclina- 
tion angles v of the orbital planes of the stars and their az- 
imuthal velocities Vaz- A similar colour table as in Fig. 2 has 
been used, but its peak value at red refers to a density which is 
a factor of about 3 lower than in Fig. 2. The dashed lines indi- 
cate the loci of stars with W = ±50 km/s and Vaz ranging from 
20 to 270 km/s as expected for thick disk stars. The locations 
of stars assigned to halo overdensities are marked by yellow 
symbols. 

Fig. 4. Three-dimensional rendering of the wavelet transform 
of the phase space distribution of nearby halo stars. The ren- 
dered surface corresponds to the outer countours in Fig. 2. 

Table 1. Parameters of the proposed halo star streams 
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1.5 
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Fig. 3 as a strong overdensity at v = 6°. A second feature can 
be identified on the left side in the bottom right panel of Fig. 2. 
Also this feature is seen very clearly as an overdensity on the 
left side of Fig. 3 at v = 170°. The three proposed streams (la- 
belled Si, S2, and S3) are described here for the first time, to 
the best of our knowledge. Moreover, we have found signs of 
Kapteyn's moving group as defined by Eggen (1996), although 
it is not very prominent. Its position in phase space is given to- 
gether with that of the other streams in Table 1. Note however 
that Kapteyn's star itself is not a member of our sample. 

3.2. Monte Carlo simulations 

In Table 1 we give the numbers of stars which we assign to 
each overdensity in phase space by our procedure. These are 
overplotted as yellow symbols in Fig. 3. The numbers of stars 
in each phase space overdensity are quite small. In order to as- 
sess the reality of the overdensities we have performed Monte- 
Carlo simulations of a smooth velocity distribution of halo 
stars. We adopted a triaxial Schwarzschild distribution with the 
velocity dispersions determined by Arifyanto & Fuchs (2005). 
The size of the simulated sample was set equal to the observed 
number of stars with non-disk kinematics. In a typical Monte - 
Carlo simulation the phase space distribution is rather smooth, 
but does show overdensities due to Poisson noise. We have 
run in total 100 Monte Carlo simulations and have calculated 
the residuals of the wavelet transforms of each individual sim- 
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Fig. 5. Three-dimensional rendering of the signal - to noise ra- 
tio of the wavelet transform of the phase space distribution of 
the halo stars at a level of S/N = 2.8, which corresponds to 
a confidence level of 99.5%. The expected Poisson noise has 
been determined with Monte Carlo simulations of a smooth 
Schwarzschild distribution (cf. section 3.2). 

ulation against the very smooth superposition of all simula- 
tions. For each cell in the Vaz - ^ae - v cube we have deter- 
mined the mean and variance of the residuals. Finally we have 
divided each value of the wavelet transform of the observed 
phase space distribution of halo stars shown in Figs. 2, 3, and 
4 by the square root of the variance of the numerically sim- 
ulated Schwarzschild distribution. This procedure gives then 
the S/N ratio of the detected overdensities. In Fig. 5 a three- 
dimensional rendering of this ratio at a level of S/N = 2.8 is 
shown. At azimuthal velocities less than Vaz < -220 km/s the 
variances of the residuals are nominally so small, that the S/N 
ratio becomes artificially high. Since Fig. 1 shows absolutely 
no clustering in this region of phase space, we have masked 
out this part in Fig. 5. Four of the star streams proposed here 
are clearly discernible, and are detected at a confidence level of 
99.5% or better. Kapteyn's moving group does not show up at 
this significance level. 

3.3. Orbital parameters of the halo streams 

We have determined the orbital parameters of the halo streams 
proposed here and list them in Table 1. Since some of the orbits 
are fairly eccentric, Dekker's (1976) approximation of Galactic 
orbits is not accurate enough to give precise estimates of the pa- 
rameters. This does not, however, diminish its power to project 
even such orbits close together into the same region of phase 
space, if the orbits are similar to each other! The mean guid- 
ing center radii and inner and outer turning radii given in Table 
1 have been calculated assuming a logarithmic Galactic grav- 
itational potential (t>(r) = V^^^ ln(r/l kpc) which corresponds 
to a flat rotation curve. The mean guiding center radius of a 
stream ro is determined by the minimum of the effective po- 
tential <l>eff('') = 'J>(i') + 0.5 ■ L^/r^ and the turning radii are 
given by the condition E - OefrCrpen, ?"apo)> where E denotes 
the orbital energy. As can be seen from Table 1 stream HB is 
at present close to its perigalacticon, whereas the other streams 
dive deep into the inner Galaxy. It would be intriguing if distant 
counterparts of the locally detected streams could be identified. 
For instance, Martmez-Delgado et al. (2007) have presented a 
numerical simulation of the Sagittarius stream which suggests 
that its leading arm passes in the vicinity of the Sun out of 
the direction of the North Galactic Pole nearly perpendicular 
through the Galactic plane, although this was questioned by 
Newberg et al. (2007). The stream labelled S3, which we pro- 
pose here, has precisely these kinematics. However, its orbital 
parameters deduced here are difficult to reconcile with the sim- 
ulated morphology of the Sagittarius stream. Obviously this de- 
pends to a large degree on the details of the model adopted for 
the Galactic potential. Observations as presented here might be 
helpful to constrain these. 
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